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Abstract 

We consider the level-crossing problem in a three-level system with non-linearly time-varying 
Hamiltonian (time-dependence t~^). We study the validity of the so-called independent crossing 
approximation in the Landau-Zener model by making comparison with results obtained numeri- 
cally in density matrix approach. We also demonstrate the failure of the so-called "nearest zero" 
approximation of the Landau-Zener level-crossing probability integral. 
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I. INTRODUCTION 



The Landau-Zener (LZ) model has been widely used for studying the dynamics of two- 
level quantum systems Ij. When the Hamiltonian of a time dependent system is diagonal, 
the energy levels would in general cross each other. In the presence of non-diagonal terms 
("perturbations") the crossing is, however, avoided leaving a time- varying non- vanishing gap 
between the levels. If the gap is sufficiently large, transition between the states is strongly 
suppressed and the system behaves adiabatically remaining in its particular energy state. 
However, if the levels get close to each other the adiabaticity may be violated and a transition 
between the states can take place in the vicinity of the points where the diagonal elements 
of the Hamiltonian (also called diabatic energies) cross. 

The exact analytical solution of the wave equations for this kind of systems can be found 
only in some special cases. The LZ model offers an approximate method for determining 
the probability (-Plz) of the transitions for two-level crossings by using a quasi-classical 



approach, essentially the WKB approximation. Originally the LZ-model was app 
transitions in diatomic molecules, but it has later found applications also in nuclear 



ied to 
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atomic [3| physics, quantum optics [4] and laser-driven atoms in time-dependent electric 



7,^. Recently the LZ-model has been 

n 

9[ and Bose-Einstein condensates [lOj ■ 



and magnetic 6|] fields, as well as in neutrino physics 
employed also in problems related to nanostructures 

In contrast to the two-level case, no general approximate solution, i.e. a counterpart of 
the LZ model, has been found for the level crossing problem of multilevel systems with three 
or more states. That is, there is no general rule which would estimately tell the state of the 
multilevel system after the system has passed the level- crossing region when the state of the 
system before entering the region is given and which would take into account the effect of 
level crossing. Such theoretically derived rules exist only for some special cases (see e.g. ll| 
and references therein). 

There are also some results which lack mathematical justification but are verified with a 
high accuracy by numerical studies 1^ . One such result concerns with a three-level system, 
where diabatic energies behave all differently as a function of time so that energy levels 
cross with one another (see Fig. 1). It has been shown numerically that the probability for 
transitions in this kind of systems can be obtained as a product of the relevant two-state LZ 
probabilities given that the Hamiltonian describing the system has a linear time dependence. 
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This is what one would expect in the case when the crossing points of diabatic energy levels 
are well separated, since then only two states at a time are close to each other while the 
third one is separated from these two with a larger energy gap. What is more surprising is 
that this simple rule is shown to be valid even when all the crossing points are close to each 
other That is, the so-called independent crossing approximation works not only in the 
case where the crossings are well separated but also in the case where they overlap. 

In this paper we will study this kind of three-level system in the case of a nonlinear 
time dependence. We have come across with such a situation when studying evolution of 
neutrino states in environments where the time dependence of the Hamiltonian is due to 



13|. The matter affects different neutrino 



a varying density of the background medium 
types differently, which results in level crossing phenomena. These effects depend on the 
density distribution of the matter, and they are typically non-linear. We will consider 
dependence, which is a typical situation in supernovae. This implies the time dependence 
as light neutrinos are highly relativistic. We will solve the evolution of the three-level 
system numerically by using the density matrix approach, which we suppose to give the 
correct description of the behavior of the system. We then make a comparison of the results 
of the numerical study with the results obtained by the LZ theory and the independent 
crossing approximation. 

II. THREE-LEVEL NEUTRINO SYSTEM 

We will consider a system of the electron neutrino (z/g), another active neutrino (i^a,^ = 
/ior r), and a so-called sterile neutrino (z/^). These neutrinos behave differently in a medium: 

encounters both charged and neutral current weak interactions, neutral current inter- 
actions only, and has no interactions with matter at all. The interactions with a medium 
appears as a "potential energy" (Ve, K, K), which is proportional to the density of the mat- 



ter p and the coupling constant of the weak force Gi? js]. The density p varies as a function 
of the location of the neutrino, or the distance r the neutrino has traveled from its creation 
point. 

The evolution of the system is described by the Schrodinger equation (the Hamiltonian 
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FIG. 1: Diabatic energy levels of a system consisting of the electron neutrino (the steepest straight 
line), a sterile neutrino (the next steepest line) and the muon or tau neutrino (the horizontal line) 
in a medium of varying density with normalization defined in the text. Non-diagonal elements of 
the Hamiltonian makes the system to avoid the crossings making the adiabatic states evolve as 
indicated by the different line formats in the figure. 



is assumed real for simplicity) 
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Hes Hss + Vs Hsa ' (l) 

Hea Hsa Haa + Va J \^a j 

The elements i^bc, which constitute the time-independent part of the Hamiltonian, corre- 
spond to the Hamiltonian in vacuum, the non-diagonal elements (6 7^ c) arising purely from 
the mixing of neutrino flavours. The neutrino mixing has been observed in experiments 
through neutrino oscillation phenomena among the neutrinos created in the atmosphere by 
cosmic rays and among neutrinos originating in the nuclear processes taking place in the 
core of the Sun. In the following we will not fix our consideration to any particular physical 
situation but analyze the three-level problem generically. 



The HbcS are the elements of the Hamiltonian in vacuum presented in the (t'e, J^s, t'a) 
flavor basis. They are obtained from 



U- 



2E 



ml 



\ 



(2) 



U 



\ 



(3) 



y ml J 

where the matrix U can be parameterized in terms of three mixing angles 6i,i = 1,2,3 as 
follows: 

/ 

C1C3 S1C3 S3 

-S1C2 — C1S2S3 C1C2 — S1S2S3 S2C3 
\ S1S2 - C1C2S3 -C1S2 - S1C2S3 C2C3 y 
where q = cos 6i and Si = sin 6i . We will consider the mixing angles and neutrino energy as 
free parameters and we fix their values as is suitable for our purposes. 

The mass terms in the diagonal elements of the Hamiltonian result from the relativistic 
approximation Ei = ^Jp"^ + mf ~ p + m1/2E once the irrelevant overall phase exp(ip) is 
omitted. In other words, the behavior of the states in the resonances depends only on the 
energy differences between the states, not on the absolute values of energies. Similarly, one 
can also normalize the potential energies. As is common, we will take out the neutral current 
contribution in a form of an overall phase, leaving the following effective matter potentials 
or the electron neutrino, the other active neutrinos, and the sterile neutrino, respectively 



Va = 



(4) 



V 



1 

72 



GpNe. 



We have assumed ordinary matter where N^. oc p, and we have for simplicity assumed that 

The starting point of the LZ approximation is an analytic continuation of the classical 
action integrals into the complex coordinate plane. In our case this results in the following 
transition probability between the two states Vb and z/c (6, c = e, s, a) (in natural units 



Pbf = exp 



-2Im J ° dt^UHbb + V) - {H,, + + 4if, 



(5) 



Here to is a point (a branch point) on the upper half-plane of the complex time, the one 
nearest to the real time axis, where the quantity under the square root vanishes. The 
integrand is the energy gap between the two eigenstates (quasi-energies) of the 2x2 sub- 
Hamiltonian spanned by the states Uf, and u^. Of course, the eigenstates of the submatrices 
generally differ from the proper energy eigenstates of the full 3x3 Hamiltonian. 

In the case where the time dependence is linear the integral in ([5]) can be easily solved 
and the probability is given by the following simple expression 

= exp[-7r7bc/2], (6) 

where the parameter is proportional to the square of the perturbation connecting the 
states b and c, i.e. of Hbc, and inversely proportional to the slope difference of the unper- 
turbed energy levels at the crossing point tr: 

In contrast to the linear case, in non-linear cases the integral ([5]) cannot in general be 
presented in closed form but one has to evaluate it numerically or approximate it with a 
suitable series expansion We will evaluate the integral numerically. 



III. RESULTS 



We now describe the results of our study of the validity of the independent crossing 
approximation in a three-neutrino case with a non-linear time dependence of energy levels. 

We first compute numerically the time evolution of neutrino states using the density 
matrix formalism. That is, we solve the matrix equation 

p = i[p,H], (8) 

where p is the density matrix. The diagonal elements pee, Pss and paa of the density matrix 
give the probability for the system with a given initial state at high positive density to end 
into the states u^, Vs and i/q, respectively, at high negative density. Negative density is, of 
course, just a formal notion; physically, a neutrino at negative density is interpreted as an 
antineutrino at positive density as the matter potential V has an opposite sign for neutrino 
and its antiparticle. (In the cases we will consider it will actually be enough to compute 
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the evolution to zero density as no level crossing will take place in negative densities). The 
relevant contribution to the level crossings is achieved in a limited range where the matter 
potentials Vf, {b = e, s, a) are of the order of the perturbations iff,c- The quantity that 
characterizes the size of this range is the transition width 6rhc, which we define as follows: 

Hbc 



5r^ 



be 



(9) 



{d/dr){V,-V,)\r=r^ 

When one is well outside this range practically no level crossing takes place but the system 
remains in the adiabatic energy state. 

The (i^e, ^s, ^a) content of the eigenstates of the Hamiltonian of course varies due to the 
mixing of the neutrinos. The varying ceases when the potential energy terms Vb overrule 
the other terms Hbc of the Hamiltonian and the system ends to one of the states Ug, Vs and 
Pa at large negative densities. As no level crossing takes place at negative densities in our 
case, the flavour states at large negative densities have one-to-one correspondence with mass 
eigenstates at zero density. 

By adjusting the values of the parameters of the Hamiltonian one can move the locations 
of the crossing points in respect to each other. We have used this in order to test the 
dependence of the independent crossing approximation on the distance of the crossing points. 
As far as the crossing points are well separated, that is, their transition width ranges are far 
from overlapping, the crossings can be clearly considered as independent two-level systems 
covered by the Landau-Zener model. Then the probability for the level crossing is given by 
the product of the Landau-Zener probabilities of the two-level crossings on the route. 

When the crossing points are so close to each other that their corresponding widths 
overlap, the situation is different and it is not a priori justified to consider the two-level 
transitions as independent. Nevertheless, numerical studies show that in the case of linear 



time-dependence, the independent crossing approximation still gives the correct result |12 |. 
We have verified the validity of this result for three-neutrino system in the case of linearly 
varying matter density. In the following we will study the non-linear case. 

We assume that our system consists in the beginning at high matter density from the 
adiabatic state z/3 = i/g purely, and we calculate the probability for having the system in the 
adiabatic state vi in the end at zero density (corresponds to the state at large negative 
density). As one can deduce from Fig. 1, the transition from the adiabatic state z/3 to 
the adiabatic state vi requires non-adiabatic behaviour of the system in two crossings it 
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encounters, i.e. in the 21- and 32-crossings. In the LZ model the probabihty P^i is thus 
given by 

Pif = PifP,f. (10) 

In Fig. 2 we present the probability of the 3 — > 1 transition calculated by the LZ ap- 
proximation {P-[i^) and numerically using the density matrix formalism {P^^i) for different 
distances between the two crossing points. We have changed the locations of the crossing 
points by varying the parameter Am^, which is the difference between the adiabatic energy 
states z/i and 1^2 at zero density (see Fig. 1). At high values of Am^ the crossing points are 
well separated and their corresponding transition regions do not overlap. With decreasing 
Am^ the crossing points come closer to each other, and at Am^ = 1 x 10~^eV^ the two 
transition regions start to overlap. At Am^ = 0.5 x 10~^eV^^ the regions are equally wide 
and overlap completely. At Am^ = 0.2 x lO^^eV"^ they are again fully separated. 

0.35 1 1 , 1 1 1 1 , 1 1 



0.3 - 




0.05 - 

qI 1 , 1 1 1 1 , 1 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

Am^[eV^] X10-' 

FIG. 2: The probability of the 3^1 transition as a function of the parameter Am^ as calculated 
by using the LZ theory and the independent crossing approximation (solid line) and numerically 
using the density matrix formalism (dashed line). The two transition regions overlap fully at 
Am2 = 0.5 X lO^^eV"^. 



8 



As the plot in Fig. 2 shows, the Landau-Zener model seems to work in our non-linear 
case relatively well, though underestimating the transition probability to some extent. Any- 
way, the overlapping of transition regions seems not to dramatically affect the validity of 
the independent crossing approximation^. The situation is actually better than the Fig. 2 
indicates, as we will now explain. 

Fig. 2 shows a curious behaviour of when Am^ is decreased beyond the maximal 
overlap region. The LZ result seems to be in a serious contradiction with the numerical result 
of the density matrix approach around the point Am^ = 0.2 x 10~^eV^. This behaviour is 
actually an artifact of simplistic use of the LZ theory and can be understood as follows. As 
mentioned, in the integral ([5]) the upper bound of the integration is a point on the upper 
half-plane of the complex time nearest to the real time axis, where the quasi-energy gap 
between the two states vanishes. In general, the closest branch point gives the dominant 
contributions and the contributions of the other branch points with a larger imaginary part 
are exponentially small, as assumed in the LZ theory. 

In our case it happens, however, that at Am^ = 0.2 x lO^^eV^ there are two branch 
points ((to)i and (^0)2) with an equal imaginary part, and actually the role of the closest 
branch point moves from one to to the other at this value of Am^. That is, different to gives 
the dominant contribution for large values of Am^ than for small values. This is illustrated 
in Fig. 3, where we plot P^^^ for both of these two zeros. Obviously, the contributions of 
the both zeros (actually, of all zeros) should be taken into account when evaluating P-^i^ . It 
goes beyond the scope of the present work to find out an appropriate way to combine the 
different contributions coherently (for a discussion on this matter for a two-level systems, 
see jl^). Nevertheless, it is quite obvious that once this summation is properly done the 
odd behaviour of the P^i^ in Fig. 2 will level off and that the the result obtained by the LZ 
theory and independent crossing approximation gets closer to the numerical result obtained 
using the density matrix approach than the Fig. 2 shows. 

^ The LZ probability (l5l) n eeds to be modified in the extreme inadiabatic hmit. We have done it according 
to the ansatz of ref. [1J| 
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FIG. 3: The contributions corresponding to the two dominating branch points (to)i ^-nd (^0)2 to 
the transition probabihty P^i^ (dash-dotted and dashed Unes, respectively). Also shown are the 
numerical result obtained by the density matrix approach (thick solid line) and the incoherent sum 
of the contributions from the two branch points (thin solid line). 

IV. CONCLUSIONS 

We have studied in this paper the vahdity of the independent crossing approximation for 
a three-level system with non-linear time dependence in the framework of the Landau-Zener 
model. In particular, we have investigated the case where the transition regions of two level 
crossings overlap, and we have found that considering the crossings as separate still gives a 
reasonably good approximation for transition probability as compared with the numerical 
estimate obtained by using density matrix approach. We have demonstrated also that the 
simplistic use of the Landau-Zener theory, where the tunneling integral is calculated only 
for the "nearest zero" in the complex time plane, leads in some occasions to very inaccurate 
results. 
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